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ABSTRACT 


A digital filter representation of the ASQ-81 
magnetometer is derived from the s-plane transfer functions 
of the system through the use of a bilinear transformation, 
A FORTRAN computer program is written which applies this 
representation to time-sampled total magnetic field data in 
order to obtain a time series representation of ASQ-81 
filtered total field, A series of simulations and a field 
experiment are conducted which verify the program output. 
Applications of this program include usage in conjunction 
with geomagnetic field data in order to produce a new data 
set representative of geomagnetic noise observed by Navy MAD 
(Magnetic Anomaly Detection) aircraft with the potential to 


investigate techniques of reducing geomagnetic noise in MAD 


aircraft. 
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I. INTRODUCTION 


The detection and location of submarines (and other 
magnetic bodies) through the discrimination of changes. or 
anomalies in the Earth's magnetic field is called Magnetic 
Anomaly Detection or MAD. In this technique, a magnetometer 
measures the magnitude of the Earth's magnetic field and 
provides an indication of that magnitude, or, more usually, 
an indication of changes in the magnitude of the Earth's 
field. These changes, or anomalies, can indicate the 
presence of magnetized bodies which may or may not be a 
submarine. 

The magnetometer currently in use in the United States 
Navy for use in this MAD process is the AN/ASQ-81 
metastable helium vapor total field magnetometer. 

Research is currently being conducted at the Naval 
Postgraduate School in Monterey, California, in various 
aspects of the applications of magnetometers, including 
Magnetic Anomaly Detection (MAD). Within the context of 
this research, magnetic field measurements are made through 
the use of sets of wire wound coils vice any specific 


magnetometer or magnetic detecting system. The data 


collected through the use of these coils is evaluated and 


ee a Fe 
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Pe processed in a variety of methods for different project 


goals. 


Ae ake Boh, 


. This thesis project is designed to produce an acceptable 
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alternative to the physical presence of an experimental 


' 


AN/ASQ-81 magnetometer at the postgraduate school by 


e 


paar 


allowing the determination, in conjunction with other 


ats 


Ai a 
wobec t ata tate 


research in progress, of the output of the AN/ASQ-81 


. 


ton Wy irre 


magnetometer from the data collected from the _ school's 


oy 

m measurement coils. It is hoped that this will assist future 
us 

“ research projects as, for example, in allowing a 
3 determination of environmental noise of such characteristics 
a as to affect the AN/ASQ-81 magnetometer operationally with 
‘ the eventual goal of providing an environmental noise index 
= ° or a system of removing such noise from the magnetometer- 
ss detection system. 


II. GEOMAGNETICS REVIEW 


A. EARTH'S MAGNETIC FIELD 
1. Constituents of the Geomagnetic Field 

The most common method of specifying the consituting 
parts of the geomagnetic field is to divide the field in 
terms of distance from the center of the Earth. This 
method results in three classifications: internal, 
crustal, and external. [Ref. 1] 

The internal field originates in the core region and 
is the most stable field, containing only extremely low 
frequency temporal variations. The crustal, or anomalous, 
field arises from modifications made on the internal field 
by materials and structures in the Earth's crust. These 
variations are not constant with regard to spatial 
locations, and comprise part of what is known as geological 
variations. The external field is the most dynamic and 
arises from many sources, including the interaction betwren 
the solar wind and the Earth's magnetic field. 

In addition to this method of defining the Earth's 
magnetic field is the method of time variations. This 
method consists of considering that part of the field which 


varies with periodicities greater than about one year as the 
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steady field and everything else as the variation field. 
[Ref. 2] 

The steady field consists of the internal field, also 
refered to as the main field. Slow variations of the main 
field with periods of years or longer are referred to as 
secular variations. 

There are various elements that contribute to the 
geomagnetic field, some of which are external to the Earth's 
surface. External contributions make up only a small part 
of the steady field, but play an important role in the 
variation field. These external sources include current 
systems in the Earth's upper atmosphere affected by solar 
electromagnetic radiation and gravitation, solar corpuscular 


radiation and the interaction of solar plasma with the main 


field, and the effect of the solar inteplanetary field. 


ety 
atte 
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[Ref. 3] 


The geomagnetic field changes with time. As 
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previously mentioned, very slow variations in the main field 
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with periods of on the order of years to thousands of years 
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are referred to as secular variations. Secular variations 
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are caused by a variation in the strength or orientation of 
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the Earth's center dipole. 


* Other time variations of the field can be categorized 
os , : baer ; : ree . 

Me iato quiet variation fields and disturbed variation fields. 
mY 

ao Disturbed variation fields include geomagnetic 
By micropulsations, which are of particular interest to 
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operational forces as these can mask target signatures and 
are therefore a source of noise to MAD sensors. 

Quiet variation fields are those which are not due to 
disturbances in the interplanetary environment and which 
vary slowly and regularly. [Ref. 3] 

Disturbed variation fields are geomagnetic field 
variations that appear to be the result of interplanetary 
environmental changes and do not posess a simple 
periodicity. These variations include ionospheric 
disturbances, the aurora, geomagnetic storms, and 
geomagnetic micropulsations. 


2. Elements of the Magnetic Field Vector 


The geomagnetic field vector is characterized at any 
point by its direction and magnitude. This is commonly 
accomplished through a system of coordinates as shown in 
Figure 2.1. The field is measured in terms of local 
coordinates with respect to true North. [Ref. 3] 


The various coordinates are referred to as magnetic 


elements and are defined as follows: 


B: Total field intensity (the symbol F is 
sometimes also used, as in this figure.) 

H: Horizontal component 

X: Northward, or NorthSouth component 

Y: fastward, or EastWest component 


Z: Downward, or Vertical component 
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Figure 2.1 Magnetic Field Elements [Ref. 4]. 


D: Declination or magnetic variation 
This is the angle between X and H and 
is measured positive eastward. 
I Inclination or Dip Angle. 
This is the angle between H and B (or F) 


and is measured positive downward. 


III. THE AN/ASQ-81 MAGNETOMETER 


A. DESCRIPTION OF SYSTEM OPERATION 

The Magnetic Anomaly Detecting set currently in use in 
the US Navy is the AN/ASQ-81 magnetometer. This set is 
used to locate and classify submerged submarines by sensing 
disturbances in the Earth's magnetic field (anomalies) 
caused by the presence of the magnetic mass of the 
submarine. The disturbance of the Earth's field is detected 
by the magnetometer, processed through filtering circuits, 
and amplified. The output signal of the magnetometer is 
displayed on a chart recorder for interpretation by an 
operator. 

The magnetic detecting set is a metastable helium vapor 
Magnetometer. The operation of the magnetometer is based on 
the light absorbtion properties of helium gas subjected to 
certain light stimulus (optical pumping), radio frequency 
excitation, and the Earth's magnetic field. The 
magnetometer consists of a helium lamp, lens and polarizer 
to generate a beam of polarized light radiation. This 
focused and polarized light beam is directed through a 
helium absorbtion cell to an infrared (IR) detector. Sone 
of the helium gas in the absorbtion cell is maintained in a 


metastable state by application of VHF excitation. 
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Figure 3.1 : Metastable Helium Magnetometer [Ref,5]. 


The Earth's magnetic field imposes a magnetic force upon 
the excited helium vapor atoms to force the atoms into one 
of three energy sublevels. This is called the Zeeman 
effect. The rate or frequency of atomic precession caused 
by this effect is called the Larmor frequency. A helium 
lamp is used to optically pump the atoms in the absorbtion 
cell, with the result that the polarized light energy 
passing through tne absorbtion cell will polarize 
(magnetize) the helium atoms in the absorbtion cells by 
selectively pumping the Zeeman levels of the energy of the 
helium atoms in the cell. The magnetization direction is 
determined by the polarization of the photons from the 


helium lamp. 


RF energy is then introduced to the absorbtion cell in 
the form of an additional magnetic field imposed through the 
use of coils oriented perpendicular to the precessed 
polarized helium atoms in the absorbtion cell and energized 
by a variable frequency RF oscillator. The RF oscillator is 
tuned to the Larmor frequency, which results in 
depolarization of the atoms. The atoms attempt to equally 
repopulate the Zeeman energy levels. However, the helium 
lamp is still beaming polarized light energy into the 
absorbtion cell, causing the atoms to absorb light energy 
and rise to an excited energy levbel. This absorbtion of 
light energy is detected through the use of an infrared 
detector. The RF oscillator frequency producing maximum 
light absorbtion is called the resonant frequency, and is 
determined through the use of a servo loop from the infrared 
detector to the RF variable frequency oscillator. 

Therefore, any change in the Earth's magnetic field 
intensity will result in a change in the Larmor frequency of 
the helium atoms in the helium absorbtion cell. This new 
Larmor frequency will be detected by the ASQ-81 
magnetometer. Since the gyromagnetic ratio of helium is 
28.024 HZ per gamma, this detection of the resonant 

on frequency provides a measurement of the Earth's magnetic 


ee field intensity at any given time. A change in the Earth's 
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p magnetic field intensity could signal the presence of a 


t submerged submarine. 
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The output resonant frequency developed by the 


“4 a. ¥ 
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magnetometer is converted to a proportional output voltage 


a, 


which is filtered through the Magnetic Anomaly Detection 
= (MAD) bandpass filters for environmental noise reduction 


cs and utilized to drive a chart recorder for observation by 


7” 
“ an operator. [Ref. 6] 

x B. TRANSFER FUNCTIONS 

- Transfer functions for the AN/ASQ-81 filters were 
E obtained from the manufacturer of the AN/ASQ-81 detecting 
3 set, Texas Instruments of Dallas, Texas. These transfer 
e ‘ functions are listed in Appendix A and are in the form of 
- H(s), that is, the frequency domain, or S domain, where 
¥ S = jw. The s-domain representation for transfer functions 
S is routinely utilized to express output system 
x characteristics for given system inputs. As the S domain 
x representation is not utilzed further in this discussion, it 
wd will not be further explained. 

RS As the output signal of the ASQ-81 magnetometer is 
< filtered through a fixed high-pass system, then through a 


selectable low pass system and a selectable high pass system 


listed in this order. 


ey (as shown in Figure 3.2 below), the transfer funtions are 
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Further discussion will be made of the selectable filters 


later. 
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IV. DIGITAL FILTERING MODELLING OF SYSTEMS 


A. SEQUENCE REPRESENTATION OF TIME FUNCTIONS 
1. Signal Representation 

A signal can be defined as a function that conveys 
information, generally about the state or behavior of a 
physical system. Although signals can be represented in 
many ways, the information conveyed by the signal is 
contained in a pattern of variations of some form. Signals 
are represented mathematically as functions of one or more 
independent variables, one of the most common of which is 
time. 

The independent variable of the mathematical 
representation of a signal may be continuous or discrete. 
Continuous time signals are signals that are defined over 
continualy values of time and are therefore represented by 


continuous-variabled functions. Discrete time signals are 


defined at discrete time intervals and are therefore 


represented by functions whose independent variable(s) take 


on discrete values only. Discrete-time signals are 
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represented as sequences of numbers. [Ref. 7] 


In addition to the fact that the independent 
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variables can be either continuous or discrete, the signal 


amplitude can be either continuous or discrete. Digital 
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oo signals are those for which both time and amplitude are 
an discrete. Analog signals are those for which both time and 
ay auplitude are continuons. 
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Digital signal processing deals with transformations 
of signals that are discrete in both time and amplitude, 
a usually represented by sequences of numbers. The nth number 


in the sequence x being processed is usually represented as 


ma x(n), and is formally written as: 

Ce 

EY x=[x(n)], -O0< n < +00 

Bats In general, an arbitrary sequence can be expressed as 

ine 00 

os x(n) = x(k) d(n-k) 

As ag k= -0O 

f ‘ 

Bs where d(n-k) is the unit sample at time k, In other words, 

wiv 

Ate an arbitrary sequence may be expressed as a sum of scaled, 

_ shifted unit samples, where the scaling factor is equal to 

os the amplitude of the sequence at that time, 

ae 2. Linear Shift-Invariant Systems 

‘Ne 

pot 

‘eau A system is defined mathematically as a unique 
transformation or operator that maps an input sequence 
{x(n)] into an output sequence [y(n)]. This is denoted 

a as: 

peat 

2 y(n) = T[x(n)] 

pre 


and is often depicted as in Figure 4.1. 
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Figure 4.1; Representation of Transformation of an Input 
Sequence to an Output Sequence. [Ref.7] 
Classes of discrete time systems are defined by 
placing constraints on the transformation T[]. 
The class of linear systems is defined by the 
principle of superposition. If y (n) and y (n) are _ the 
1 2 
responses when x (n) and x (n) are the inputs, then a system 
1 2 
is linear if 
Tlax (n) + bx (n)] = aT[x (n)] + bDI[x (n)] 
1 z 1 2 
= ay (n) + by (n) 
1 ps 
for any arbitrary constants a and b. This, together with 
the concept of representing a sequence by a sum of delayed 
and scaled unit-sample sequences, suggests that a linear 
system can be characterized by its unit-sample response. 
Specifically, let h (n) be the response of the system to 


k 
d(n-k), a unit sample occurring at n=k. Then 


co 


y(n) = T[ x(k) d(n-k)] or, 


= -00 


22 


oe ae ae TS a A Lt a 


y(n) = x(k) T[d(n-k)] = x(k) h (n) 
k= -0O k= -@O ; 
Thus the system response can be expressed in terms of the 
response of the system to d(n-k). 

The class of shift invariant systems is characterized 
by the property that if y(n) is the response to x(n), then 
y(n-k) is the response to x(n-k), where k is a positive or 
negative integer. When the index n is associated with time, 
shift-invariance corresponds’ to time-invariance. The 
property of shift invariance implies that if h(n) is the 
response to d(n), then the response to d(n-k) is simply 


h(n-k). Therefore oo 


y(n) = x(k) h(n-k) 
k= -00 
and any linear shift-invariant system is completely 
characterized by its un’ :-sample response h(n). 
A subclass of .ineir shift-invariant systems consists 
of those systems for which the input x(n) and the output 
y(n) satisfy an Nth-order linear constant-coefficient 


difference equation of the form 


N M 

3 a y(n-k) = b x(n-r) 
Canned k r 

k=0 r=0 
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If the assumption of causality is made about’ the 
system, a linear difference equation provides an explicit 
relationship between the input to the system and the output 
of the system. This can be seen by rewriting the previous 


equation as 


N M 
y(n) = ) c y(n-k) + » d x(n-r) 
k r 
k=] r=0 
where c =-a /a andc =b /a. 
k k 0 r r 0O 
Thus the nth value of the output can be computed from 
the nth value of the input and the N and M past values of 
the output and input, respectively. The difference equation 
not only represents the system for theoretical purposes, 
but it may also serve as a computational realization of the 


system. The z-Transform makes use of this property to 


realize systems. 


B. THE z-TRANSFORM 
l. Description of the z-Transform 
The z-transform plays an important role in the 


analysis and representation of discrete-time linear shift- 


invariant systems. The z-transform, X(z), of a sequence 
x(n) is defined as CO 
-n 
X(z) = x(n)z 
n= -®&© 
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where z is a complex variable. This representation of the 
z-transform is referred to as the two-sided z transform. The 
one sided z-transform consists of the same summation for 
terms of n greater than or equal to zero. For the case 
that x(n)=O for n<0, the one sided and t'» sided z 
transforms are equivalent. 

By expressing the complex variable z in polar form as 
zs Ber the z-transform can be interpreted as the Fourier 


transform of x(n) multiplied by an exponential sequence. 


For r = 1, that is, for Izl = 1, the z-transform is equal to 


2a 


A eee, 


the Fourier tranform of the sequence, 


wor 


. 2. The Bilinear Transformation 


The transfer functions of analog systems are most 


. ae 


often expressed in terms of s = jw (see section III B.). 
This corresponds to the analog frequency response of the 
system. This analog frequency response can be "mapped", 
that is, tranformed to the z-plane from the s-plane through 
the use of the bilinear transformation. The eff<cct of 
utilizing the bilinear transformation is to convert a system 
transfer function in terms of the variable S into the system 


transfer function in terms of the variable z, The 


transformation itself is: 


a 
(2/T) + 

os Ui: Seeae eS 
7. (2/T) - s 

s 
a where T is the sampling period, that is, the time between 
2 data samples. 

oe Thus a transform can be made from one plane to’ the 
ES other. In this manner, the transfer function, H(z), of a 
s system may be obtained. 
x The bilinear tranformation equations may be shown to 
= hold in general, and the use of this transformation may be 
= shown to yield stable digital filters from stable analog 
% filters [Ref. 7]. The bilinear transformation maps. the 
on 
ce . imaginary jw axis in the s-plane onto a unit circie (of the 
i region of convergence) in the z-plane, with the left half 
- s-plane mapped onto the region inside the circle and _ the 
= right hand (region of instability) s-plane mapped onto the 
sc region outside this circle [Ref. 8]. A complete discussion 
i of the z-transform is available in several texts, some of 
~ which are listed in the Bibliography. 

C. THE DIGITAL COMPUTATIONAL ALGORITHM 
% In implementing a digital filter on a digital computer 
such as the IBM 3033, the input-output relationship of the 
— signals through the system being synthesized must be 
= converted to a computational algorithn. The algorithm is 
Ee specified in terms of a set of basic computations of 
mn elements. For the implementation of discrete-time systems 
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described by linear constant coefficient difference 
equations, such as the AN/ASQ-81, it is convenient to choose 
as these elements the basic operations of addition, delay, 
and multiplication by a constant. The computational 
algorithm for implementing the filter is then defined by a 


structure or network consisting of an interconnection of 


these basic operations. For a system transfer function of 
the form M 
~k 
bz 
k 
k=0 ¥(z) 
H(z) Ss £4ss ee eoewZoeveo = —— 
X(z) 
1- a z 
k 
k=1 


the difference equation relating input and output is easily 
written down directly from the system function and is given 

by N M 
y(n) = >: a y(n -—- k) + > b x(n - k) [Ref. 7] 

k k 

k=1 k=0 
This difference equation can be interpreted directly as a 
computational algorithm in which the delayed values of the 
input are multiplied by the coefficients b , the delayed 
values of the output are multiplied by ie coefficients 
a , and the resulting products are added. It is now easy 
to see the process to be followed in obtaining the 


computational algorithm for the AN/ASQ-81 magnetometer 
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transfer function. The z-transform of the system transfer 


function is obtained through the use of the bilinear 
transformation, and is then converted into a difference 
equation relating input and output signals, thence to a 
FORTRAN computer program. A table of z-transforms of system 
functions is included in Appendix B. 

In the FORTRAN computer program realization of the tctal 
system computational algorithn, each filter block is 
transformed into a separate difference equation and 
algorithm. This was done to enable a "building block" type 
approach to the program, and to minimize computational and 


roundoff errors. 


D. THE CASCADE FORM OF THE COMPUTATIONAL ALGORITHM 

Even though the direct form realization of the digital 
filter design may be perfectly satisfactory in a theoretical 
sense, it may be less than desirable in the context of 
realization through the use of a general purpose computer of 
fixed register length. The parameters of a digital filter 
are usually obtained with a high degree of accuracy, which 
results in a faithful realization of the desired system. 
When these parameters are quantized, as in a finite memory 
register within a computer, the frequency response of the 
resulting digital filter may differ appreciably from the 


original design. In fact, the quantized filter may fail to 
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meet design specifications although the unquantized filter 
does. [Ref. 7] 

The sensitivity of the filter response to errors in the 
filter parameters is dependent upon the structure of the 
filter realization. Therefore, in the event of an 
unacceptable change in the frequency response of the filter 
due to quantization errors, it is often possible to minimize 
the effect of these errors through an alternate filter 
realization structure. An alternate structure to the 
previoulsy discussed direct form realization is the cascade 
form realization. 

The direct form network structures were obtained directly 
from the system function H(z) written in the form of a ratio 
of sums. If this ratio is factored into a product of 


polynomiais of the form 
[(N+1)/2] “1 % 


1k 2k 
H(z) = A | fo wee een ene ee 
aa =3 
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this product represents a general distribution of poles and 
zeros and suggests a set of structures consisting of a 
cascade of first and second-order subsystems. There is 
considerable freedom in the design of the subsystems, but it 
is best to realize the systems using a minimum of storage. 
The expression of H(z) in this form indicates the 


presence of poles and zeros in pairs. If poles and zeros 
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are not present in pairs, one of the coefficients B or a 
2k 2k 
will be zero as appropriate. An implementation of such a 
cascade structure with the use of minimum memory can be 
obtained through a direct form II realization of each second 
order subsystem using techniques similar to the direct form 
implementaticn utilized previousiy. A cascade realization 


of a sixth-order system, such as the ASQ-81 system, using a 


direct form II realization of each second order subsystem 


would appear as in Figure 4.2 below. 


Cascade Structure With a Direct Form II 
Realization of Each Second Order Subsystem. 
[Ref. 7] 


Figure 4.2 


There is, theoretically, considerable flexibility in the 
Manner in which the poles and zeros are paired together and 
in the order in which the resulting second-order subsystems 
are cascaded. However, although all such pairings and 
orderings are equivalent for infinite~precision arithmetic, 
they may differ considerably in practice owing to finite 


word length effects of roundoff and truncation. 
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Figure 5.1: Output of First Stage Filter of Digital Filter 
Computer Program With Sinusoidal Input in Simulation, 
Unfortunately, the second stage output of the filter 

showed an instability within the program design, indicated 

by the output of the filter being a sinusoid of increasing 


magnitude, as indicated in Figure 5.2 below. 
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Figure 5.2: Output of Second Stage Filter Design With Input 
of a Sinusoid. 
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The stability of the third stage of the filter design 
6 was investigated by inputting the sinusoid directly to the 
third filter, and found to be stable. A check of the 
derivation of the equations, coefficients, and programming 
steps of the second (unstable) filter of the design failed 
to indicate the cause of the instability. 


Computation of the poles of the z transfer function, 


ie H(z), of the second stage of the filter confirmed the 
Re instability of the design. The poles were computed to be: 
= 0.92 + 0.12181, 1.07 + 0.13401, 0.8611, and 1.1564. Of 
ad these six poles, three lie outside the region of convergence 
BS for the z-plane, that is, within the unit circle discussed 
at 

= previously in Chapter IV. 

e : The second stage of the filter was therefore 
s redesigned using the cascade form of the direct form 
i realization (direct form II), and tested in simulation. A 
im copy of the software used in the simulation is enclosed in 
a Appendix F, 

a The output of all three filter stages of the program 
ss were stable, as indicated in Figures 5.3 through 5.7 below. 
z The amplitude decrease and phase shift expected were 
8 observed. The "damped overshoot" of the second stage output 
a4 is due to the fact that, for values of the input function 
: prior to time zero in the simulation, utilized in the input- 
x output signal difference equations for the filter, the input 
we Signal was ‘set at QO. This resulted in an instantaneous 
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change of the input signal from 0 to 


introduced in the simulation at time O+. 
the filter is the filter's attempt 


instantaneous jump in magnitude of the 


the finite value 
The "overshoot" of 
to "match" this 


input signal. When 


the input signal to the filter in the simulation is zero 


at time zero, this overshoot effect does not occur. 
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Figure 5.3: Input Signal to Digital 
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Sinusoid of Frequency 0.1 HZ and Amplitude + 1. 
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Figure 5.4: Output of First Stage of Digital Filter Program. 
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: Pigure 5.5: Output of Second Stage of Digital Filter 
* Progran. 
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Figure 5.6: Output of Third Stage of Digital Filter Program. 


This simulation was run with inputs of sinusoids of 
various frequencies in order to check the stability of the 
; filter design at frequencies throughout the operating range 
of the AN/ASQ-81 magnetometer. In all cases, the design was 
stable, and the expected amplitude changes and phase shifts 

occured. 

2. Noiselike Inputs 

The simulation was also run with inputs consisting of 


a sinusoid of a frequency which should be passed through the 
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AN/ASQ-81 added to sinusoids of frequencies which should 


4 


have been filtered by the magnetometer and random noise. The 


filter performed as expected, with the sinusoid of a 


ta ee, 


a 
2 passable frequency passed by the filter, and spurious noise 
: and sinusoids attenuated severely. The results of a 
= 
4 
f simulation consisting of a sinusoid of passable frequency, a 
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filterable sinusoid, and uniformly distributed random noise, 


all of amplitude + 1, are presented in Figures 5.7 through 


5.10 below. 
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Figure 5.7: Input to Filter ~ 0.1 HZ Sinusoid, 10 HZ 
Sinusoid, Uniformly Distributed Random Noise of Amplitude 
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Figure 5.8: Output of First Filter Stage. 
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Figure 5.10: Final Filter Stage Output. 


As can be seen, the digital filter program succeeds 
in filtering out random noise and signals of frequency 
components above the band pass of the magnetometer. 

In order to ensure that the digital filter 


representation of the magnetometer has the same amplitude 
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versus frequency characteristics of the AN/ASQ-81 


Magnetometer, a simulation program was written which inputs 
sinusoids of varying frequencies and computes the Root Mean 
Square (RMS) value of the filter output and the signal 
input, then computes the decibel (dB) attenuation of the 
filter at that frequency. A copy of this program is 
included in Appendix G. A plot was made of the dB 
attenuation versus frequency for the filter and compared 
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Figure 5.11: Plot of Attenuation Versus Frequency for 
Sinusoidal Inputs for Digital Filter and Analog Filter. 
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with the measured frequency performance of the AN/ASQ-81 
Magnetometer, which was supplied by Texas Instruments, Inc., 
and does not include the effects of the fixed high pass 
filter. Consequently, the data shown in Figure 5.11 is a 
comparison of the data supplied by Texas Instruments and the 
output of the test program, which also does not include the 
fixed high pass filter. As can be seen, the performance of 
the filter is extremely similar to that of the magnetometer 


itself. 


3. Anderson Function Simulations 

The next step in the simulation phase was’ the 
introduction to the filter of Anderson function simulations. 
The shape of the signal amplitude of the output of a 
magnetometer passing through the sphere of influence of a 
magnetic anomaly (submarine) is a function of the dip angle 
of the geomagnetic field, the magnetic heading of the track 
of the magnetometer (or the aircraft), the magnetic 
heading of the anomaly (submarine) dipole, and the 
lateral range between the magn:tometer (aircraft) and the 
anomaly, Anderson functions [Ref. 9] are mathematical 
representations of three basic components of signals which, 
when taken in various linear combinations, describe the 


shape of these anomaly signals. The equations for the 


Anderson functions are: 


2 5/2 
(First Anderson Function) f = I1/(1 +B ) 


w oO 
it 


where (velocity)x(time) 

range at CPA 
or, a dimensionless parameter defined as the distance 
traveled along the magnetome er (aircraft) track divided by 
the slant range at closestpoint of approach (CPA) 
Bx f 
1 0 2 


Bxf =B xf 
Z 1 0 


(Second Anderson Function) f 


(Third Anderson Function) f 


The Anderson functions were introduced into the 
filter program in a noise-free signal environment in order 
to observe the output signal and ensure that it was a "MAD- 
like" signal. A rigorous determination of the actual output 
signal would have been extremely difficult to obtain, so a 
comparison was made with the output of a computer simulation 
program provided to NPS by Mr. Joe Rice of Texas 
Instruments. When the sampling rate of the program was 
adjustec to equal that of the Texas [nstruments program, 8 
HZ, the two program outputs were observed to be very 
Similar. The Anderson function simulation inputs and 
outputs of the program are depicted in Figures 5.12 through 
og. Bis The Texas Instruments program outputs were obtained 
in the form of time series plots of discontinuous data 
points, and were therefore not conducive to replotting for 


comparison. 
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Figure 5.14; Second Anderson Function Input. CPA at Time 
2.5 Minutes. 
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Figure 5.15: Filter Output for Second Anderson Function. 
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Figure 5.16; Third Anderson Function Input. CPA at Time 
2.5 Minutes. 
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Figure 5.17: Filter Output for Third Anderson Function, 


The filter output for all three Anderson function 
inputs did appear to he "MAD-like" signals, and did closely 
resemble the Simulation output obtained from Texas 
Instruments, Inc. 


4, Impulse Function Response 
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" The response of the filter program was also observed 
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+. when the input was a unit impulse function. Again, the 
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output was compared to that of the Texas Instruments' 
computer program. The outputs of the two programs were 
observed to be, again, very similar, as can be seen in 
Figure 5.18, where the response of the NPGS filter is 
represented by a solid line and that of the Texas 


Instruments filter by a chain-dash line. The abrupt "jumps" 
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of the Texas Instruments response are due to the translation 
of the output plot supplied to this plot. The plot supplied 
by Texas Instruments was, again, discontinuous points of 


poor resolution, and it was necessary to interpolate values 


45 


co OD ge Ce ee a ora ne Re Neti Tes ee SOME PMN TRC EL 8 RAM Be Ue my M ! nn Ste lee a he eR ae eae eee) at eee ating SE | Sa ge hag hs ae 


Ee a a er Re peg! Me OO ome a eg Neg OR Bare NL RR eR RR ae aS we Ciena ee ae 


in order to generate Figure 5.18. This resulted in the 
broken appearance of the plot. Even so, the similarity of 


the outputs can be observed. 


B. EQUIPMENT SETUP 

Following the simulation phase of the experiment, actual 
magnetic field measurements were introduced to the filter in 
order to test the response of the filter. Magnetic field 
measurements were made at the La Mesa field test site near 
the Naval Postgraduate School in Monterey. The output of an 
AN/ASQ-81 magnetometer, a Schonstedt magnetic field sensor, 
and the school's coil sensor, oriented along the Earth's 
magnetic field, were pulse code modulated (PCM) and 
transmitted via VHF radio to recording devices at the 
Postgraduate school. The recording of a two hour long data 
collection period was transferred to digital data tape for 
use by the school's IBM3033 general purpose mainframe 
computer. 

In the first test of the digital filter program, the 
output of the Schonstedt sensor, which represents 
fluctuations of the Earth's total field, was used as_ the 


input to the computer program. A comparison of the output 


of the computer program, with this approximation to the 


total field fluctuations as input, to the output of the 


wy rg 


Pasa a ea a tee ee Oe 
Feta ORME fe ete fe ee OE 


AN/ASQ-81 should provide an indication of the proper 


functioning of the computer filter program. The results of 


we eee 


the test are shown in Figures 5.19 through 5.21 on the pages 
following. Figure 5.19, the Schonstedt sensor output, shows 
several instances of PCM dropouts, that is, occasions where 
the pulse code modulation signal was not correctly read by 
the computer for some reason. At such occurances, the data 
point value used by the computer is a random number and does 
not reflect the true value of the data. The problem with 
these PCM dropouts is that the computer does not recognize 
them as invalid data points and will use them in 
computations. This can (and does) cause problems in the 
computation of Fourier transforms, spectral characteristics, 
etc. Additionally, this will also impact the _ proper 
functioning of the digital filter program which is the 
subject of this thesis. PCM dropouts are visible at times 
6, 8, 142, and 220 through 226 seconds on the plot of the 
Schonstedt sensor output. An examination of Figure 5.20, 
the filter program output, reveals the programs attempt to 
"follow" these PCM dropouts. It should be recalled that 


previous simulations indicated the filter's tendency to 


2 "follow" sudden changes in the input signal, with a 
% relaxation time required for the filter to steady out. This 
= effect is apparent in the output of the filter program at 
a times corresponding to those of the PCM dropouts in the 
a Schonstedt sensor’s time series plot. It can be seen that 
Be 

: : this overshoot tendency resulted in an output significantly 
= different from the actual AN/ASQ-81 output at these times. 
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a 
= If the PCM dropout induced differences are neglected, it can 
a be seen that the shape of the output of the filter program 
Re is remarkably similar to that of the AN/ASQ-81 magnetometer, 
me : although noisier. Note that the output of the AN/ASQ-81 
ay magnetometer exceded the maximum voltage amplitude which the 
= pre-amplifiers of the data collection system were able to 
handle and resulted in a truncated signal from time 40 to 
Ns time 60 seconds. It can still be seen, however, that the 
= filter program output is very similar to the time signal 
te which would have been displayed without this truncation. 
“4 
Be It should be noted that the amplitudes of the time series 
eo signals of the progrem output and the AN/ASQ-81 magnetometer 
= differ consideradly. In the case of the program output, the 
x peak amplitudes are on the order of 1.4 nanoteslas along the 
a vertical scale, while the peak amplitudes of the output of 
= the AN/ASQ-81 magnetometer are on the order of 0.7 
z nanoteslas along the vertical scale. This is because the 


bay. input signal to the filter program is an approximation to 
the total field difference time series signal, and some 
amplitude difference could reasonably be expected. The 


intent of this initial test was to investigate the output 


time series shape, and an exact correlation was not 
oe! ‘i expected. It is worth noting that the digital filter 
ow 
we program will perform its function on any time series signal, 
regardless of units. This means that a signal may be 
51 
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operated upon either before or after conversion from 
whatever units it was originally measured to magnetic field 
strength units. 

Therefore it appears that the digital filter program is 
functioning properly. When a close approximation to. the 
fluctuations of the total field time series signal is used 
as the input to the computer program, the output of the 
program is similar to the time series output of an AN/ASQ~81 
magnetometer. 

The final stage in the testing process was a conversion 
of the time series output voltage signal of the coil antenna 
sensor, which was aligned along the Earth's magnetic field, 
into a total field fluctuation time series representation 
for the same time period as before, and then to use this as 
the input to the digital filter progran. A comparison of 
the resultant time series output of the program with the 
actual AN/ASQ-81 magnetometer output would validate the 
proper functioning of the program. 

Conversion of the time series antenna sensor output 
voltage signal into total field fluctuations in nanoteslas 
was accomplished through the use of a computer program 
designed by Capt. Kurt Stevens, USAF, a student at the Naval 
Postgraduate School, as his Master's thesis [Ref 10]. The 
output voltage time series is stored in an array, then a 
Fourier transform is performed on the stored data, resulting 


in the Fourier spectrum of the data. This spectrum is 
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corrected for the characteristics of the coil antenna sensor 
to obtain the Fourier spectrum of the total field data. A 
reverse Fourier transform gives the time series signal for 
total magnetic field in nanoteslas. 

This time series signal was used as the input to the 
digital filter program and compared with the output of the 
AN/ASQ-81 magnetometer. Figures 5.22 through 5.24 show the 
raw coil antenna data, the total field time series data, and 
the program output time series for a 6 minute period of the 
test. Figure 5,22 shows the raw coil antenna data series. 
The number of PCM dropouts should be noted, as these will 
influence the performance of the filter program. Figure 
5.23 shows the computed total field time series. Note that 
the PCM dropouts evident on the raw time series plot are 
evident on the computed total ficl- time series plot also, 
and thus inputted to the filte program as valid data 
points. Additionally, there are two "jumps" in the plot of 
total field fluctuation (Figure 5.23) which are also 
inputted to the filter program as valid data points. These 
"jumps" are located at 128 and 256 seconds and are caused by 
the method of processing blocks of data for the conversion 
to total field fluctuation. A block of 128 seconds of data 
is processed at a time, and the results of each block are 
scored in an array. This results in a slight amplitude 


difference between the last data point of one block and the 
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first data point of the next block of data. This slight 
difference is manifested as a signal jump. 

A comparison of Figures 5.24 (program output) and 5.21 
(AN/ASQ-81 output) show that the filter program gives a time 
series output very similar to that of the actual 
magnetometer. The first 20 seconds of the program output is 
somewhat dissimilar to that of the AN/ASQ-81, due either to 
the initial “start up" delay of the filter program or to 
distortion of the total field fluctuation time series, 
There is a PCM dropout at time 11 seconds which contributed 
to the distortion, 

Following this, however, it can be seen that the program 
output is very similar to that of the magnetometer, except 
at 128 and 256 seconds, which show the effects of the false 
signal jumgs caused by the total field fluctuation 
conversion. There is also a noticeable four to five second 
time delay between the AN/ASQ-81 output and that of the 
filter program. As this delay is not evident in a 
comparisen of the AN/ASQ-81 output and that cf the filter 
program with the Schonstedt sensor as the input, it can be 
inferred that this time delay is caused either by the 
program which converts the raw coil data to total field 
fluctuation data, or by a phase (and hence time) change of 
the voltage signal due to the coil sensor itself. A 
comparison of the raw coil data in Figure 5.22 to the 


converted coil data in Figure 5.23 indicates no time shift, 
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and hence the deduction can be made that there is a time 
delay inherent within the coil sensor itself. 

Other than the differences «© the four to five second 
time delay and the distortions caused by the false signal 
jumps, the program output is extremely similar to the output 


of the AN/ASQ-81 magnetometer, 
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VI. CONCLUSIONS 

The intent of this thesis was to design and test a 
digital filter computer program which would, when given a 
time series input of fluctuations in the total magnetic 
field, deliver an output time series representation of the 
output of an AN/ASQ-—81 magnetometer. This purpose has been 
realized. 

The computer program contained in Appendix I has’ been 
proven to output a time series signal which is very similar 
to that of the magnetometer. The major limitations of the 
output signal are a finite time delay of about five seconds 
between the AN/ASQ-81 magnetometer signal and the output 

. signal of the program, a sensitivity of the program to false 
data points such as those caused by PCM dropouts and false 
signal jumps caused by processing large data blocks, and the 
inherent limitations of the program caused by its dependence 
on the use of digital data tapes and the IBM 3033 mainframe 
computer. 

The five second time delay is not considered to be an 
important limitation to the program, as it was intended as a 
research tool for programs currently in progress at the 
Naval Postgraduate School. Instances where this time delay 
might become important would be in areas of simultaneous 
comparison of the program output signal with an actual 


Magnetometer, in target location algorithms using time 
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delays, or in correlation studies between different sensors. 
In correlation studies using coil sensors, the effects of 
the time delay would cancel out, as all coil outputs would 
be similarly delayed. In target location algorithms, the 
target location errors due to the time delay could be 
adjusted for simply while in computer simulation, and flight 
testing could not reasonably be accomplished without the use 
of an actual magnetometer as the’ sensor. Lastly, in a 
comparison of the program output with an actual sensor, the 
time delay can, again, be compensated for. In short, these 
limitations are not considered excessive, especially as the 
apparent cause for the delay is not the filter program. 

In the primarily intended purpose of the filter program, 
magnetic noise studies, the time delay is not considered to 
be a problem. 

The problem of false data points caused by PCM dropouts 
and signal jumps due to conversion to total field 
fluctuations is more serious. False data points cause 
inaccuracies in the output time series and could adversely 
affect later projects. Unfortunately the PCM dropout 
problem is one which is endemic to the data collection 
system presently being used at the postgraduate school, and 
not to the filter program itself. It is imperative that 
users of this program are aware of the PCM dropout problem 


and of the effects it may entail upon their specific 


Se a eee ee 


™ research, A large number of PCM dropouts in a time series 
could render that series unusable. Similarly the case of 
the false data jumps caused by conversion to total field 
fluctuations is not within the filter program. Further 
investigation of this problem is necessary in order to 
eliminate it. 

The last problem, that of reliance upon the digital data 
tape/IBM 3033 computer system, is, like the PCM problem, one 
which is not endimic to the filter program but rather to the 
data collection system being used. A change of data 
collection system may, at some future time, remove the 


~ reliance upon the PCM/digital tape/IBM 3033 data system (and 


Lee 


a a 


ae hence too the data block conversion problem which results in 
falso data jumps), but this is unlikely at this time. Users 
should be aware of this dependence and of possible effects 


upon specific research projects. 
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APPENDIX A 
AN/ASQ-81 FILTER TRANSFER FUNCTIONS 


Fixed High Pass Transfer Function: 
2 
80 S$ 
H(S) = ----5----+-----+-+----- 
80 3° + 20 S + 1 


Selectable High Pass Transfer Functions: 


A. 0.04 HZ H(S) = 
7s 2 2 
40.82834 § 45,28317 § 
eo aan on ss om ee we ee ee en X wee ee ae ee nn en ee ee 
40.82834 S$ *+ 12.52096 $ +1  45.28317 §- + 11.00999 § + 1 
2 
57.576688 § 
Xo ewe ee ee ee ee me ee eee 
57.576688 S* + 7.41498 S$ + 1 
B 0.06 AZ H(S) = 
2 2 
18.14591 § 20.12587 § 
18.14591 & 4 8.34727 S +1 20.12587 8 4 7.33999 $3.1 
2 
25.58964 § 
An SSeS SS= De mene ae Rg ye 
25,.58964 S° + 4.94332 S$ + 1 
C. 0.08 HZ H(S) = 
ae 2 2 
10.20708 § 11.32080 § 
ee ee eee ee ee ee Xo wwe we ee ee a ee ne 
10.20708 $*+ 6.26045 $ + 1 11.32080 s* + 5.50500 § + 1 
2 
14.39417 8 
x ee ee ee 


14,39417 S* + 3.70749 $ +1 
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0.10 HZ 


Zz 2 
6.53253 § 7.24531 § 
em ee eee ee gSeaA nese See = x a a ono irae mee i ea a eerie 
6.53253 S + 5.00836 S +1 7.24531 S$ + 4.40400 S + 1 
2 
9.21227 § 
Xo ee wee ee ee we ee ee ee ee 
9.21227 $* + 2.96599 S$ +1 
Selectable Low Pass Transfer Functions 
A. 0.2 HZ H(S) = 
1 | 
te ote Gt Gee ee SO he > he ee ee ee ee ee ee ee Oe ee ee ee es ee ee ee ee x Ses ery a seed BI Ae a Cea 
0.3143 32 + 1.0741 S$ + 1 0.2501 S“~ + 0.6209 S$ + 1 
B. 0.4 HZ H(S) = 
1 1 
Soe Ar Saenger eae ee 
9.07858 S* + 0.53706 S +1 0.06252 S“ + 0.31044 S + 1 
C. 0.6 HZ H(S) = 
1 1 


0.03492 S~ + 0.35804 S + 1 0.02779 g? + C.20696 & + 1 
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APPENDIX B 


AN/ASQ-81 Z TRANSFORM FILTER TRANSFER FUNCTIONS FOR DIRECT 
FORM I REALIZATION 


For fixed high pass filter: 


-1 aD 
BFHPO + BFHP1*Z + BFHP2*Z 
HZ) area teee eee sprrt etna _y-------- 


1 - AFHP1*Z.~ - AFHP2*Z 
where BFHPO, PFHP1, BFHP2, AFHP1, . AFHP2 are constants 
tabulated in Appendix D. 
For selectable high pass filter: 
BSHPO + BSHP1*Z = + anes + sages 
+BSHP4*Z~* + BSHP5*Z> + BSHP6*z ~© 
H(Z) 8 ecewnercnennz Speer Shes ese cs a en a= =: isa 
1 - ASHP1*Z ~ - ASHP2*Z ~ - ASHP3*Z 
- ASHP4*Z~* - ASHPS*z ~> - asHP6*z~° 


where, fer low frequency cutoff of 0.04 HZ: 


BSHPO = C.99471378 

BSHP1 = -5.9682827 ASHP] = 5.9894021 
BSHP2 = 14.920707 ASHP2 = -14,947051 
BSHP3 = -19.894276 ASHP3 = 19.894225 
BSHP4 = 14,920707 ASHP4 = ~14,894327 
BSHP5 = -5,9682817 ASHPS = 5.9472141 
BSHP6 = 0.99471372 ASHP6 = -0.98945296 
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For Low frequency cutoff of 0.06 HZ: 


BSHPO = 0.9920813 


BSHP1 = -5.9524928 ASHP1 = 5. 9841070 
BSHP2 = 14.881232 ASHP2 = -14.920650 
BSHP3 = -19.841643 ASHP3 = 19.841528 
BSHP4 = 14,881232 ASHP4 = -14.841757 
BSHPS = -5,9524927 ASHPS = 5.9209919 
BSHP6 = 0.99208212 ASHP6 = -0.98422128 
For low frequency cutoff of 0.08 HZ: 
BSHPO = 0.98945806 
BSHP1 = -5.9367483 ASHP1 = 5.9788144 
BSHP2 = 14.841871 ASHP2 = -14.894276 
BSHP3 = -19,789161 ASHP3 = 19.788959 
BSHP4 = 14,841871 ASHP4 = -14.789365 
BSHPS = -5.9367476 ASHPS = 5,8948841 
BSHPE : 0.98945802 ASHP6 = -0.97901720 
: For low frequency cutoff of 0.1 HZ: 
o BSHPO = 0,98684156 
= BSHP1 = -5,9210493 ASHP1 = 5.9735244 
us BSHP2 = 14,802623 ASHP2 = -14.867941 
BS BSHP3 = -19.736831 ASHP3 = 19.736516 
Re BSHP4 = 14.802623 ASHP4 = -14.737149 
is BSHPS = -5.9210491 ASHP5 = 5,8688889 
i BSHP6 = 0.9868415 ASHP6 = -0.97384065 
: 
= 
= 
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BSLFO + BSLP1*Z + BSLP2*Z + BSLP2*Z + BSLP4*Z 
H(Z) = ------------- fo een ee te ee ee a a dos 
1 - ASLP1*Z - ASLP2*Z - ASLP3*Z - ASLP4*Z 


-2 


where BSLPC, BSLP1, BSLP2, BSLP3, BSLP4, ASLP1, ASLP2, 


ASLP3, ASLP4 are constants tabulated in Appendix D. 
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APPENDIX C 


AN/ASQ-81 Z TRANSFORM FILTER TRANSFER FUNCTIONS DIRECT 
FORM SI REALIZATION 


For fixed high pass filter; 
ee We -1 ~2 
BFHPO + BFHP1*Z + BFHP2#*Z 
UC) ia ee eee ese ese es cee 
1 - AFHP1*Z71 - aFHP2*Z~2 


where BFHPO, BFHP1, BFHP2, AFHP1, AFHP2 are constants 


tabulated in Appendix D. 
For selectable high pass filter: 
-1 


H(Z) = ASHPL x -~-------~---4-----------5 x 
1 - ASHP3*Z~~ - ASHP2*Z 


-1 -2 -1 -2 
1- 2*#2 +2 1 - 2*2Z + Z 
Oe ee ce ec eae Ye ae = OS eae TG aN an VE Ce ROTA =2 
1 — ASHP4*Z ~ ASHP5*Z 1 ~ ASHP6*Z - ASHP7*Z 
I where ASHP1, ASHP2, ASHP3, ASHP4, ASHP5, ASHP6, ASHP7 are 
hs 
MS constants and tabulated in Appendix D. 
ed For selectable low pass filter: 
os Or, oe ay) ~3 -4 
ee BSLPO + BSLP1*Z + BSLP2*Z + BSLP3*Z + BSLP4*Z 
H(Z) 3 -nenn nnn nr enn ee ee yen en ne ee ee eee 
3 


oa eS 


where BSLPO, BSLP1, BSLP2, BSLP3, BSLP4, ASLP1, ASLP2, 


OF ES Eh a as 
Re Srna 
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ASLP3, ASLP4 are constants tabulated in Appendix D. 
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APPENDIX D 
Z TRANSFORM REALIZATION DIFFERENCE EQUATIONS 


With reference to Figures 3.2 and 4.2, the following 
difference equations are used to model the AN/ASQ-81 
magnetometer filter transfer functions. The input to the 
fixed high pass filter is called SIG(I), where I is’ the 
current data sample. The output of the fixed high pass 
filter, which is the input to the selectable high pass 
filter, is YO(I), and the output of the selectable high 
pass fiiter, the input to the selectable low pass filter, 
is called YPO(I). The output of the filter is called 
ASQ(I). (I-1) denotes a time delay of cne sample, and so 
forth, and the symbol * denotes multiplication. 

For the fixed high pass filter: 

YO(I)=BFHPO*SIG(I) + BFHPI*SIG(I-1) + BFHP2*SIG(I-2) + 
AFHP1*YO(I-1) + AFHP2*YO(I-2) 

where: 

BFHPO = 0.9980499222938581 


BFHP1 = ~1.9960998445877161 AFHP 1 1.9960983216843922 


BFHP2 = 0.998049922238581 AFHP2 = -0.9961015674910398 
For the selectable high pass filters: 

XI(I) = ASHPI*YO(I) + ASHP2*XI(I-2) + ASHP3*XI(I-3) 
XII(I) = XI(I) + XI(1-2) - 2*xXI(I-1) 

XIITI(I) = XIICI) + ASHP4*XITII(I-1) + ASHP5S*XIII(I-2) 


XIV(I) = XIIICI) -~ 2*XIIICI-1) + XIII(I-2) 


XV(I) = XIV(I) + ASHP6*XV(I-1) + ASHP7*XV(I-2) 

YPO(I) = XV(I) - 2*®XV(I-1) + XV(I-2) 

For the low frequeacy cutoff at 0.04 HZ: 

ASHP] = 0.9947137/89347288 ASHP2 = -0.9952196910157882 
ASHP3 = 1.9952137256322473 ASHP4 = 1,9962028201103847 
ASHP5 = -0.9962082013015601 ASHP6 = 1.9979855321466768 
ASHP7 


-0.9979897681491607 


For the low frequency curo £ at 0.06 HZ; 


ASHP] = 0.99208212771993253  ASHP2 = -0.9928381306174365 
ASHP> => 1.9928247245317958  ASHP4 = 1.9943056083414792 
ASHP5 = -0.9943177045263504 ASuP6 = 1.9969766455640039 
ASHP7 = -0.9969861717656565 

For the low frequency cutoff at 0.08 HZ: 

ASHP1 = 0.9894580558875787 ASHP2 = ~0.9904622611337722 
ASHP3 = 1.9904384565912725  ASHF4 = 1.9924092959984783 
ASHPS = -~0.9924307799269113 ASHP6 = 1.9959666590811389 
ASHP7 = -0.9959835866201267 

For the low frequency cutoff at 0.10 4Z: 

ASHP] = 0.9869415560096681  ASHP2 = -0.9880929619779491 
ASHP3 = 1.9880549117849344 AfHP4 


1.9905139074397893 
ASHP5 = -0.9905474442556225 ASRP6 


1.9949555824611562 
ASHP7 = -0.9949820174650785 
For the selectable low pass filters: 


ASQ(T) 


tt 


ASLP1*ASQ@(I-1) + ASLP2*ASQ(I-2) + ASLP3*ASQ(I-3) 


+ 


ASLP4*ASQ(I-4) + BSLPO*YPO(I) + BSLP1*YFO(I-1) 
+ BSLP2*YPO(I-2) + BSLP3*YPO(I-3) + BSLP4*YPO(I-4) 
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For the 


BSLPO = 
BSLP1 = 
BSLP2 = 
BSLP3 = 
BSLP4 


For the 
BSLPO = 
BSLP1 = 
BSLP2 = 
BSLP3 = 
RSLP4 = 
For the 
BSLPO = 
BSLP1 = 
BSLP2 = 


BSLP3 
BSLP4 = 


high frequency cutoff at 0.2 


0.0000000452615229 
0.0000001810464917 
0.0000002715697375 
0.0000001810464917 
0.0000000452616229 


high frequency cutoff at 0.4 


0.0000006918001209 
0.0000027672004837 
0.0000041508007256 
0.0000027672004837 
0.0000006918001209 


high frequency cutoff at 0.6 


0.0000033463317975 
0.0000133853271900 
0.000020077990785G6 
0.0000133853271900 
0,0000033463317975 


Sgt tt ate te 


ASLP1 
LSLP2 
ASLP3 
ASLP4 


ASLP1 
ASLP2 
ASLP3 
ASLP4 


ASLP1 
ASLP2 
ASLP3 
ASLP4 
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Ba cs) gs fa 


se 


3.9082436339591027 
~5.7285022156249328 
3.7321935213310065 
-0.9119356638511430 


HZ: 


3.8173771378993420 
-5.4670046844062743 

3.4812554457127576 
~0,.8316389680077603 
HZ: 


3.7274299052305002 
-5.2152772583906819 
3.2462216641520829 
~0.7584278523006615 
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DIGITAL SOFTWARE FOR SIMULATION - DIRECT FORM (SINUSOIDS AS INPUT) 
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DIGITAL SOFTWARE FOR COMPUTATION OF SYSTEM AMPLITUDE VERSUS FREQUENCY 
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